New limits on variation of the fine-structure constant using atomic dysprosium 
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We report on the spectroscopy of radio-frequency transitions between nearly-degenerate, opposite- 
parity excited states in atomic dysprosium (Dy). Theoretical calculations predict that these states 
are very sensitive to variation of the fine-structure constant, a, owing to large relativistic corrections 
of opposite sign for the opposite-parity levels. The near degeneracy reduces the relative precision 
necessary to place constraints on variation of a competitive with results obtained from the best 
atomic clocks in the world. Additionally, the existence of several abundant isotopes of Dy allows 
isotopic comparisons that suppress common-mode systematic errors. The frequencies of the 754- 
MHz transition in ^^*Dy and 235-MHz transition in ^®^Dy were measured over the span of two years. 
Linear variation of a is found to be a/a = (— 5.8±6.9) x lO"^'^ yr~^, consistent with zero. The same 
data are used to constrain the dimensionless parameter fee , characterizing a possible coupling of a 
to a changing gravitational potential. We find that ka — (—5.5 ± 5.2) x 10~^, essentially consistent 
with zero and the best constraint to date. 

PACS numbers: 06.20.Jr, 32.30. Bv 



Variation of fundamental constants was first formu- 
lated by Dirac as the Large Numbers hypothesis [U [5] . 
The observation that dimensionless ratios of quantities 
such as the age of the universe to atomic time scales 
and the electromagnetic to gravitational force between 
a proton and electron were of the same order of magni- 
tude, ^ 10"*", led to the hypothesis that these ratios were 
functions of the age of the Universe. A consequence of 
this hypothesis is a gravitational constant, G, that scales 
inversely proportional to the age of the universe. Al- 
though modern experiments based on lunar ranging [3| 
have ruled out present-day variation of such magnitude, 
the variability of fundamental constants remains an ac- 
tive area of theoretical and experiment research. Any 
such variation would be a violation of the Einstein Equiv- 
alence Principle (EEP) and an indication of physics be- 
yond General Relativity (GR) and the Standard Model 
(SM) of particle physics [11|5]. 

Changing constants would manifest in a wide range of 
physical observables. The dimensionless electromagnetic- 
coupling constant, the fine-structure constant, a, is of 
particular importance due to the implications of its vari- 
ation on atomic clocks and time keeping. Any variation 
of a would lead to a change in the relative frequencies of 
co-located clocks even in the absence of external fields. 
This is forbidden by an assumption of EEP. In this letter 
we report new constraints on variation of a with respect 
to time and changing gravitational potential from a com- 
parison of radio-frequency transitions in two isotopes of 
atomic dysprosium (Dy) [SI [7]. These new results im- 
prove on our earlier constraints [El [9] by almost two or- 
ders of magnitude and are competitive with existing lim- 
its from other experiments pUHlT) . 



The most stringent laboratory constraints on variation 
of fundamental constants come from clock-comparison 
experiments. We restrict our attention to clocks based 
on transitions in atoms and molecules. The ratio of any 
two such clock frequencies can be written as |18j 
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where A ]s a, dimensionless factor dependent on atomic 
structure, /ie = "T-e/"^p is the electron-proton mass ra- 
tio, and /iq — 'niql Kqcd is ratio of the quark mass to 
QCD-mass scale. The dimensionless constants ^e and /i^ 
are important for comparisons involving transitions with 
hyperfine structure [101 I16j or molecular transitions |19) . 
The sensitivity coefficients, Ka,e,q, depend on the par- 
ticular frequency ratio under consideration. A summary 
of coefficients for various comparisons can be found in 
Table m 

In Dy we make use of an 'accidental' degeneracy of 
energy levels to greatly relax the measurement precision 
necessary to place competitive limits on variation of a. 
Large relativistic corrections to electron energies in Dy 
create an almost complete degeneracy of opposite-parity 
excited states, labeled A and B by convention (Fig. II]). 
This system has been the subject of investigations span- 
ning over two decades, including an attempt to measure 
parity nonconservation [5D1 [H] . Recently an analysis of 
the data from the present work has also been used to 
place stringent limits on violations of Lorentz symmetry 
and the Einstein Equivalence Principle |22] . 

The energy difference corresponding, i/^a = {^b — 
tj^jh, is sensitive only to variation of a. In practice, 
however, any measurement must have a standard 'ruler' 



A Even 






J = 9 



Odd 



\ 1 .4 nm 



669 lira 



4.8nmV J = 9 

\; 4f%s6p 



J = i 



564 nm ,' 



4f5d6s= 



833 nm 



/ 

"■"Dy 




/ 

/ 754 MHz 

j 






\ 235 MHz 


'«Dy 

A: J =10 
4f''5d6s 

\ T ~ 8 us 


B: J= 10 

4f5d-6s 
T > 200 us 



J=8 

4f°6s^ 




FIG. 1: Partial energy diagram for Dy showing states of 
interest. Preparation of atoms is accomplished via two laser 
excitations and a spontaneous decay with 30% branching ratio 
into metastable state B. Atoms are excited from state B to 
A by a resonant, frequency-modulated rf electric field. State 
A decays with lifetime ~ 8 /xs. A photomultipler tube and 
lock-in amplifier detect the 564-nm fluorescence. The bottom 
right inset shows typical lock-in signals for ^^^Dy at the 1st 
and 2nd harmonics of the modulation frequency. 



for comparison. The frequencies in our experiment are 
measured with respect to the stabilized oscillator of a ce- 
sium (Cs) beam standard, which introduces sensitivity 
to variation of both /ie and /iq. Changes in the frequency 
ratio Vba/i^cs can be written 



Aln^ = if„^ 
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As shown in Table m the Dy/Cs frequency comparison 
is over six orders of magnitude more sensitive to varia- 
tion of a than to variation of /ig and /i^. This is a rela- 
tive enhancement of sensitivity, rather than an absolute 
enhancement, owing to the near degeneracy of levels A 
and B . At our present level of measurement precision, 
variation of /Xg or /.tg would only be observable at lev- 
els orders of magnitude larger than stringent constraints 
placed by other experiments |10j . Thus our experiment 
is effectively sensitive only to variation of a, 
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Instability of the Cs reference, a > 30 yr old HP5061A, 
presents another source of concern for measurements 
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TABLE I: Sensitivity coefficients for several clock compar- 
isons. CSO refers to crystal-sapphire oscillator. The large 
sensitivity of the Dy transition frequency to variation of a 
is a relative enhancement due to the near degeneracy of the 
electronic states involved in the transition. Column refer- 
ences are for experimental details. Calculations of sensitivity 
coefficients can be found in Refs. i23. 24 . 



spanning several years. A separate comparison between 
the Cs reference and a GPS stabilized Rb oscillator (Sym- 
metricom TS2700) is performed during all data collection 
as a check against this. The fractional instability of the 
Cs reference, as compared to Rb reference, has been mea- 
sured to be < 10~^^ yr~^7 well below our dominant mea- 
surement errors. The influence of Cs-reference instability 
is ignored, and the magnitude of the frequency \vba\ is 
assumed to vary with a according to [53] 



Al 



±(2 X 10^'' Hz) Aa/a, 
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where the sign is negative for Vba > and positive for 
Vba < 0. The present work is based on measurements of 
the Vba ~ 753.5 MHz and Vba ~ —234.7 MHz transitions 
in ^^^Dy and ^^^Dy (see Fig.fTj). Comparing isotopes with 
sensitivities of opposite sign allows for the cancellation of 
common systematic errors that might otherwise mimic 
variation of a in a single isotope. 

The spectroscopy is performed on a thermal beam of 
Dy atoms, produced in an oven heated to '^ 1400 K 
inside a vacuum chamber with residual gas pressure of 
^ 1Q~^ torr. After two coUimators/conductance chokes 
the atoms enter the interaction chamber where the resid- 
ual gas pressure is ~ 10~^ torr. The atoms undergo 
laser excitations at 833 nm and 669 nm, employing an 
adiabatic-passage technique [3S], followed by a sponta- 
neous decay at 1.4 /xm with 30% branching ratio to state 
B. Narrow-band lasers provide high-fidelity isotope se- 
lection. Upon excitation to state B atoms then enter 
the interaction region, where excitation from B to A 
occurs via a frequency-modulated electric field. Atoms 
spontaneously decay from state A via two steps to the 
ground state. Fluorescence at 564 nm is directed by a 
polished-aluminum light-collection system (~ 4% over- 
all efficiency) into a glass pipe, detected by a photomul- 
tiplier tube (PMT), and sent to a lock-in amplifier for 
processing. Figure [2] shows a simplified diagram of the 
experiment. 
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TABLE II: Current levels of known systematics. The total 
systematic uncertainty is obtained by adding in quadrature. 
The corresponding uncertainties for jti/aj assume two mea- 
surements separated by one year. 



FIG. 2: Schematic of the experimental set-up. Argon- 
ion lasers pump a dye laser producing 669-nm light and a 
Ti:sapphire laser producing 833-nm light. Components in 
vacuum are within the dashed boundaries, a) Skimmers col- 
limate the atomic beam, and double as conductance chokes 
for differential pumping between the oven chamber and inter- 
action chamber, b) In-vacuum linear polarizers are the last 
optical element for the laser light before interacting with Dy 
atoms, c) Lenses diverge the laser light to match the atomic 
beam divergence, d) Polished aluminum mirrors guide fluo- 
rescence to a photomultiplier tube, e) An interference filter 
wih 564-nm peak transmission suppresses stray laser and oven 
light, f) A glass pipe guides fluorescent light to a PMT for 
detection. 



The apparatus has been designed to minize the system- 
atic uncertainties presented in Table lllj In our previous 
result [8] the sensitivity was limited by the collisional per- 
tubation of energy levels by background gases j26] , poor 
suppression of Zeeman shifts owing to imperfections in 
laser-light polarization, and systematic shifts related to 
inhomogeneity of the rf field. In the new apparatus, the 
high-vacuum system reduces collisional shifts to below 
the 1 mHz level. The electric-field region has been de- 
signed to ensure field homogeneity across the range of 
operating frequencies. Doppler shifts are suppressed by 
creating an rf-standing wave in the interaction region in 
addition to orienting the k-vector of any residual trav- 
eling wave perpendicular to the atomic-beam propaga- 
tion axis. Two layers of magnetic shielding limit back- 
ground magnetic fields to below 500 /iG in all directions 
and three-axis magnetic field coils allow residual fields to 
be canceled out. 

The dominant systematic is an electronic offset in the 
acquisition electronics, which may create a shift in the 
zero-crossing of the first-harmonic signal and apparent 
shift of the transition frequency. Sensitivity to elec- 



tronic offsets is amplified by the relatively large transit- 
broadened linewidth of the transition, 7 ~ 27r x 40 kHz. 
We measure these offsets by varying the PMT gain, in 
order to change the useful signal size while leaving elec- 
tronic noise unchanged. This idea is based on the off- 
set compensation scheme presented in Ref. [27], but cur- 
rently only constrains electronic offsets at the level of 500 
nV. The importance of this effect depends on the absolute 
signal size and is reflected in the range of uncertainties 
in Table nil 

The ac-Stark shift in a two-level system is approxi- 
mately zero for a resonant electric field, with a negligible 
contribution expected from what is known as the Bloch- 
Siegert shift |7j . Strongly coupled off-resonant levels may 
lead to large shifts correlated with rf-power. A measure- 
ment of the off-resonant contributions to the dynamic po- 
larizabilities in ^^^Dy and ^^^Dy found 8v ~ 70 E"^ mHz, 
where E"^ is mean-squared field value. Typical values of 
E"^ are 4.5 (V/cm)^, corresponding to a stability of 3 
mHz/% change in rf power. The uncertainty associated 
with this systematic is conservatively estimated from an 
assumption of 10% control over the rf power in the inter- 
action region. 

Additional Stark related systematics are the dc-Stark 
effect and blackbody radiation (BBR) induced Stark 
shifts [25]. Charged particles in the atomic beam can 
cause charge accumulation on the electric field plates 
and produce DC fields. An electrode biased at 500 V is 
used to sweep charged particles out of the atomic beam, 
and the DC field is periodically measured via Zeeman- 
crossing spectroscopy [5] |5D]. These are consistently 
found to be at the level of 10 mV/cm. The temperature 
dependence of the transition frequencies has been mea- 
sured near room temperature to be -1-29(4) mHz/K and 
-34(4) mHz/K for i^^^Dy and ^^^^y, respectively The 
isotopic dependence of the sign is consistent with BBR 
induced Stark shifts, but the attribution of these shifts to 
BBR is preliminary |29j . Currently, the 2 K temperature 



stability of the interaction region is used to estimate the 
systematic uncertainty due to this effect. 

Suppression of systematics related to Zeeman shifts is 
accomplished by performing spectroscopy with the Zee- 
man structure unresolved. Linear polarizers are located 
in vacuum and are the last optical elements for the 833- 
nm and 669-mn laser light, ensuring symmetric popula- 
tion of the ±M magnetic sublevels of state B. A mag- 
netic field then leads to a broadening of the unresolved 
line, but no shift. A measured residual Zeeman shift of 
~ 2.5 Hz/mG represents a suppression of ~ 1000 from 
the sensitivity of the m = 10 sublevel. The magnetic field 
stability along the quantization axis, chosen to coincide 
with the rf field, is at the level of 20 /iG. We note that the 
magnetic field insensitive ttIb = — > rn^ ~ transition 
is forbidden between levels A and B where A J = 0. 

Residual amplitude modulation refers to a power 
imbalance of the carrier sidebands in the frequency- 
modulated spectrum of the electric field. Such an im- 
balance distorts the atomic lineshape and creates an ap- 
parent frequency shift. Poor impedance matching and 
termination of the rf transmission line made this a domi- 
nant systematic in early data at the 1 Hz level. Mea- 
suring the transition frequency with the in-phase and 
quadrature channels of the lock-in amplifier allows the 
size and stability of RAM to be measured directly with 
the atoms |5D|. This protocol was implemented begin- 
ning May 2011. In August 2011, custom narrow-band 
radio-frequency circulators (DPV CO) were acquired to 
suppress transmission-line etalons, reducing RAM to the 
level of '^ 10 ppm. The frequency shift introduced by 
this modification was measured and a correction applied 
to earlier data. 

The transition frequencies j/i64 and uiq2 measured over 
the span of two years are shown in Fig. |3] The reduced 
uncertainties beginning in May 2011 are primarily due to 
the characterization and eventual suppression of RAM. 
To constrain a linear variation of a in time, a global 
linear least-squares fit is performed, in which the two 
isotopes' data are fit by independent offsets and equal 
magnitude slopes of opposite sign. The best-fit slope of 
—0.12 ± 0.14 Hz/yr corresponds to the result 



a/a ^ (-5.8 ± 6.9) x 10"^^ yr 



(5) 



which is consistent with zero within 1 standard deviation. 
This result approaches within a factor of 3 the level ob- 
tained with the best optical clocks in the world [13] , and 
is limited by systematics. The contribution of statistical 
uncertainties is at the level of a/a ^ 1.7 x 10 
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The data are also fit by equal slopes of the same sign, 
which is sensitive to common mode systematics, but not 
variation of a. This fit gives a slope of 0.41 ± 0.14 Hz/yr. 
The 3-sigma, non-zero drift at the level of ~ 0.5 Hz could 
be explained by a drifting electronic offset which, as a 
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FIG. 3: Changes in the transition frequencies for ^^'^Dy (filled 
circles) and ^^''Dy (empty circles) over the span of two years. 
The frequencies for ^^^Dy and ^^*Dy are displayed with re- 
spect to 234,661,102 Hz and 753,513,708 Hz, respectively, a) 
The data are fit by linear functions with equal magnitude 
slopes of opposite sign (solid) and same sign (dashed), b) 
The data are fit by cosine functions with equal amplitudes but 
180° phase difi^erence (solid) and 0° phase difiterence (dashed). 
The variation of the dimensionless gravitational potential, 
scaled in relative units by 5 x 10^", is shown by the light 
solid line. 



technical systematic, is expected to be the same sign for 
both isotopes. 

Our data can also be used to constrain violations of 
local position invariance, assuming a model where fun- 
damental constants are influenced by light scalar fields 
that scale linearly with changes in the local gravitational 
potential [31]. We can express this as Aa/a — kaAU/c^, 
where AU/c^ is a change in the dimensionless gravita- 
tional potential. The ellipticity of the Earth's orbit pro- 
vides semi-annual changes in the laboratory gravitational 
potential, AU /(? — ±1.65 x 10~^°, at the aphelion and 
perihelion of Earth's orbit for plus and minus signs, re- 
spectively. To constrain fc^ the data are again fit by 
global linear least-squares to cosine functions with equal 
amplitudes but 180° phase difference. The period is equal 
to one solar year and zero phase is fixed at Earth's perihe- 
lion on Jan. 3, 2010. The best-fit amplitude of oscillation 
is found to be 0.18±0.17 Hz, providing the best constraint 
to date on the dimensionless coupling parameter [10] 



fCn 



(-5.5 ±5.2) X 10" 



(6) 



which is also consistent with zero at ~1 standard devi- 
ation. The sensitivity is again limited by systematics. 
The statistical contribution to the uncertainty is at the 
level of ka ~ 1.2 x 10~^. A global fit to the two isotopes' 
data with 0° phase difference, sensitive to common mode 



systematics, has an amplitude of — 0.17± 0.17 Hz. The 
data and best fits are shown in Figure [3] 

We have presented updated constraints on variation 
of a that represent almost two orders of magnitude im- 
provement over previous results, with the present level 
of sensitivity still limited by systematic effects. While 
more stringent control of these systematics, particularly 
electronic offsets, presents a clear avenue to achieving 
the ultimate practical statistical limit of 10~^^ calculated 
in [7], recent astrophysical evidence [35] for spatial vari- 
ation of a suggests an observable variation of a in the 
laboratory at the level of 10~^^ (331 • A new generation of 
experiments based on the spectroscopy of optical nuclear 
transitions [34j or optical transitions in highly-charged 
ions [3S] will be necessary to observe this effect. 
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